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We present a model-independent description of the exclusive rare decays B — > K*e + e~ in the 
low recoil region (large lepton invariant mass q 2 ~ m\). In this region the long-distance effects 
from quark loops can be computed with the help of an operator product expansion in 1/Q, with 
Q = {mi,, \/~cf}- Nonperturbative effects up to and including terms suppressed by A/Q and m\/m\ 
relative to the short-distance amplitude can be included in a model-independent way. Based on 
these results, we propose an improved method for determining the CKM matrix element |V^&| from 
a combination of rare and semileptonic B and D decays near the zero recoil point. The residual 
theoretical uncertainty from long distance effects in this \V u b\ determination comes from terms in 
the OPE of order a s (Q)A/mb, a 2 s (Q), m\/m\ and duality violations, and is estimated to be below 
10%. 

I. INTRODUCTION 

Radiative B decays are important sources of information about the weak couplings of heavy quarks. Experiments 
at the B factories have measured precisely the branching ratios of the exclusive rare radiative b — > S7 and semileptonic 
b — ► uev decays, and decay spectra are beginning to be probed. In addition to offering ways of extracting the CKM 
matrix elements V u b and V t d, these processes hold good promise for the detection of new physics effects (see e.g. 0). 

In contrast to the inclusive heavy hadron decays which can be reliably described using the heavy mass expansion, 
the corresponding heavy-light exclusive decays are comparatively less well understood. The theoretical ignorance of 
the strong interaction effects in these decays is parameterized in terms of unknown heavy to light B — > M form 
factors. Although lattice 0] and QCD sum rules Q have made significant progress in computing these form factors, 
they are still beset with large errors and limitations. 

In the low recoil region, heavy quark symmetry has been used to relate some of the B —> M form factors 0,0]. 
In Refs. |fj, [3 we showed that the leading corrections to these symmetry relations when nib 7^ 00 do not involve any 
non-local contributions, that is, they are characterized solely in terms of matrix elements of local operators. Here we 
show that the cancellations of non-local terms, which appear as a remarkable accident in the heavy quark effective 
theory, are easily understood by deriving the form factors relations directly from QCD at finite nib- 

For the case of b — > se + e~ decays there is an additional source of theoretical uncertainty due to long distance effects 
involving the weak nonleptonic Hamiltonian and the quarks' electromagnetic current. In B — > K*e + e~ , these effects 
are numerically significant for a dilepton invariant mass close to the cc resonance region q 2 = (p e + + p e - ) 2 ~ 10 GeV 2 . 
Usually these effects are computed using the parton model [slllOj. or vector meson dominance, by assuming saturation 
with a few lo w lying resonances ip n and using the factorization approximation for the nonleptonic decay amplitudes 
B — > K*ijj n pl [ill E2] • Such a procedure is necessarily model dependent, and its effect on the \ V u b\ determination has 
been estimated at ~ 10%. Although in principle the validity of the approximations made can be tested aposteriori 
by measuring other predicted observables, such as the shape of the q 2 spectrum or angular distributions, it is clearly 
desirable to have a more reliable computation of these effects. 

The object of this paper is to show that, near the zero recoil point q 2 ~ <7 2 lax = (tub — rriK*) 2 , these long distance 
contributions to B — ► K*e + e~ can be computed as a short-distance effect using simultaneous heavy quark and operator 
product expansions in 1/Q, with Q = {nib, y^Q 2 }- We use this expansion to develop a power counting scheme for 
the long-distance amplitude, and classify the various contributions in terms of matrix elements of operators. The 
leading term in the expansion is calculated in terms of the form factors that were necessary to parametrize the local, 
leading contribution to the decay amplitude. Moreover, the first correction, of order A/Q, is given in terms of the 
same operators introduced in Ref. [3 to parameterize the leading order corrections to the heavy quark symmetry 
relations between form factors, and is suppressed further by a factor of a s (nib). The largest second order correction, 
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of order z = m 2 /m 2 , is also calculable in terms of the leading form factors. Hence, our method for computing the long 
distance contributions introduces no new model dependencies to good accuracy. The terms we neglect are suppressed 
by m%lm\ and A 2 /m 2 relative to the short-distance amplitude, and are expected to introduce an uncertainty in \ V u b\ 
of about 1-2%. 

A model-independent determination of \V u b\ has been proposed using semileptonic and rare B and D decays in the 
low recoil kinematic region 0, 0, 0, 0] . This method uses heavy quark symmetry to relate the semileptonic and 
rare radiative B form factors. More specifically, this method requires the rare and semileptonic modes B — » A'*e+e~, 
B — > pev. D — > K*ev and D — > pev. The main observation is that, neglecting the long distance contribution to the 
radiative decay, the double ratio [T(B — > K*e + e~)/T(B — > peiA\J[T(D — > K*ev)/T(D — > pev)] is calculable since it is 
protected by both heavy quark and ST/(3)-flavor symmetries [l5(- We extend this result to include the long distance 
contributions which, as explained above, are calculable in terms of the same form factors in the endpoint region. 

The modes required for this determination are beginning to be probed experimentally. The branching ratios of 
the rare decays B -> K^H+i' have been measured by both the BABAR and BELLE [l^ (with I = e,fi) 
collaborations 

H(R _ K*f+f-) - J (°- 88 ±o.29 ± °' 10 ) x 10_6 (BABAR) 
1 ' \ (11.51^ ±0.8 ±0.2) x 10- 7 (BELLE) 

and 

B(R -> Kf+f~) - ! (°- 65 -o:i3 ± 0-04) x 10- 6 (BABAR) 
B(U->Kt I i- | ( 4 . 8 +i.O ±03±01)xl0 -7 (BELLE). 

This suggests that a determination of \V u b\ using these decays might become feasible in a not too distant future. 

The paper is organized as follows. In Sec. II we construct the operator product expansion (OPE) formalism for the 
long-distance contribution to exclusive B — > K*e + e~ decay in the low recoil region q 2 ~ <7 2 lax . This is formulated as 
an expansion in l/Q, with Q — {mb, y/q 2 }. The coefficients of the operators in the OPE are determined by matching 
at the scale Q, which is discussed in some detail in Sec. III. In Sec. IV we present the evaluation of the hadronic 
matrix elements of the operators appearing in the OPE, and explicit results for the \V u b\ determination are presented 
in Sec. V. An Appendix contains a simplified derivation of the improved form factor symmetry relations at low recoil. 

II. OPERATOR PRODUCT EXPANSION 



Woe = ^V tb V* V d{n)Qi{p) , (1) 



The effective Hamiltonian mediating the rare decays b — > se + e is 

c io 

eS = --^v tb v t * 1 £Ci(riQ 

* i=l 

where the operators Qi can be chosen as 

Qi = {s a Cf}) v -A{cf}b a ) v ^ A (2) 

Q2 = (sc) V - A (cb) V -A 

Q3 = (sb)v-A^2(qq)v-A 
1 

Qa = (Sgbp)y-A y^(g/3ffa)y-A 

1 

Q 5 = (sb) V -A^2(qq)v+A 
1 

Q& = (sgbp)y^A y^(g/?ga)y+A 



Q7 = — m b s a a^(l + 75)^^ 



g 

e 

J_ 

8tt 

e 2 



2 m b s a a^(l + IsjT^bpG^ 



Q9 = —(sb)v-A(ee) v 



8tt 2 
8tt 2 



Qw = —(sb)v-A(ee) A . 
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We denoted here (qq)v±A — 97^(1 i 75)9. The contributions of the operators Q7.9,io are factorizable and can be 
directly expressed through form factors, while the remaining operators Qis contribute through nonlocal matrix 
elements with the quarks' electromagnetic coupling jg m — J2 q QqTi^q as 

A(B -> K*e + e~) = ^V tb V*^ {(e^e)^ + {erf^e)A^} . (3) 

The two hadronic amplitudes are given explicitly by 

A ( P = -C 7 (ri^(K*(k )V )\Ma„ v q^l+ l5 )b\B(v)) (4) 

1 6 

+C 9 (n)(K*(k, 77)1^(1 - l5 )b\B(v)) - 8vr 2 - £ a(M)^ (< V, fx) 

Q i=i 

A^ = CM(K*(k,r,)\s 7fl (l-~f 5 )b\B(v)). (5) 
where we introduced the nonlocal matrix element parameterizing the long-distance amplitude 

77V) = * / d^^K^k^TQm.f^.^Biv)) . (6) 



The conservation of the electromagnetic current implies in the usual way the Ward identity (see e.g. f° r the 

long-distance amplitude 

<f77V) = 0. (7) 

Our problem is to compute T^q 2 ) in the low recoil region, corresponding to q 2 ~ m 2 . Consider the amplitude 
T~i(q 2 ) as a function of the complex variable q 2 . This is an analytic function everywhere in the complex q 2 plane, 
except for poles and cuts corresponding to states with the quantum numbers of the photon J PC = 1 . The region 
kincmatically accessible in B — > K*e + e~ is the segment on the real axis q 2 = [0, g 2 lax = (mg — my) 2 ]. 

This is very similar to e + e~ — » hadrons, which is related by unitarity to the correlator of two electromagnetic 
currents W v (q 2 ) = U(q 2 )(q^q u - q 2 g^) = i J d 4 xe lq - x (0\Tj»(0), j"(x)\Q). For this case, it is well known that at large 
time-like q 2 , both the dispersive and imaginary parts of the correlator Tl(q 2 ) can be computed in perturbation theory. 
This is the statement of local duality 20], which is expected to hold up to power corrections in l/(3[2l|,|22|. In contrast 
to e + e~ — > hadrons, the external states appearing in the definition of %(q 2 ) are strongly interacting. For this reason, 
a closer analogy is to the computation of the inclusive semileptonic width of B hadrons using the OPE and heavy 
quark expansion [23). 

The zero recoil point in B — > K*e + e~ corresponds to a dilepton invariant mass f?„ lax = (ms — mx*) 2 = 19.2 
GeV 2 and is sufficiently far away from the threshold of the resonance region connected with cc states q 2 ~ 10 GeV 2 . 
Therefore duality can be expected to work reasonably well. There arc, in addition, effects from thresholds of other 
J PC — 1 states, like the p and the T. These effects are smaller because they either enter through the operators 
Q3-Q6, which have small Wilson coefficients, or through Qf = (s q m^)(m / 36 q ) and Q% = (su)(ub) through CKM 
suppressed loops ~ V u bV* s . The effects of light states, like the p-meson, are under better control since the associated 
resonance regions are even lower than for cc. Heavier states, like the T, lie above These too are under better 

control since duality sets in much faster from below resonance than from above, as evidenced by empirical observation, 
as in the example of e + e~ — * hadrons. 

In analogy with the OPE for the inclusive B decays, we propose to expand the amplitudes %(q 2 ) in an operator 
product expansion in the large scale Q = {mb, ^fq 2 } 

v\f) - E J2 c S^/<^(°f hl (^ ( 8 ) 

fc>-2 j 

where the contribution of the operator scales like 1/Q k - The operators appearing on the right-hand side are 
constructed using the HQET bottom quark field h v , and they can contain explicit factors of the velocity v and the 
dilepton momentum q. Their matrix elements must satisfy the Ward identity Eq. Q for all possible external states, 
which has therefore to be satisfied at operator level. In addition, they must transform in the same way as under 
the chiral 5C/l(3) x SUr(3) group, up to factors of the light quark masses which can flip chirality. 

Our analysis will be valid in the small recoil region, where the light meson kinetic energy is small Ey — my ~ A. 
Expressed in terms of the dilepton invariant mass q 2 this translates into the range (mg — my) 2 — q 2 < 2tobA. In the 
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FIG. 1: Contributions to the B — > K*£ + £~ amplitude near the zero recoil point coming from different operators in the OPE 
Eq. J8J. In (a) the circled cross denotes one of the operators cK -1 >°) of the form qTiD^hy or qTgG^hy, and in (b) it denotes 
one of the 4-quark operators (qh v )(cc). The contributions in (a) are suppressed relative to the short-distance amplitude by 
A/Q (for e> (_1) ), A 2 /Q 2 (for O (0) ), and those in (b) by m 4 JQ 4 . 



particular case of B — > K*e + e this region extends about 5 GeV 2 below the maximal value g^ ax = (mB-m^.) 2 = 19.2 
GeV 2 . 

Each term in the OPE Eq. (JHJ must have mass dimension 5. The leading contributions come from operators whose 
matrix elements scale like Q 
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0[- 2) = s L {q 2 ^~q»q]h vL (9) 
£>2 = imbSLa^q u h v R . (10) 

Another allowed operator (q 2 v^ — q^v ■ q){si,h V R) can be shown in fact to scale like QA after using Eq. (|llfl . and 

is included below as Og 1 ' (see Eq. <|14f) ). These operators are written in terms of chiral fields qL.R — PL,R.q, with 
Pl.r = |(1t 75)- In the chiral limit only sl can appear, and the right-handed field sr requires an explicit factor of 
m s . 

In general, the dilepton momentum q a can be rewritten as a constant part plus a total derivative acting on the 
current 

q a (sTh v ) = (m b v a + id a )(sTh v ) , (11) 

where the two terms on the right-hand side scale like mb and A, respectively. For this reason, using q a in the definition 
of the operators gives them a non- homogeneous scaling in 1 /nib- This is not a problem in the power counting scheme 
adopted here, which counts nib and Q as being comparable. We will keep q a explicit in the leading operators Eq. JJJJ, 
(|1(J|| . which we would like to write in a form as close as possible to the short-distance operators. On the other hand, 
we expand in 1/rrib in the sub- leading operators below, and keep only the leading term in Eq. I|ll|) . 

Next we include operators whose matrix elements scale like QA. They are dimension-4 operators of the form 
qTiD^h v . A complete set of operators which satisfies the condition (J7J) and which do not vanish by the equations of 
motion can be chosen as 

(12) 
(13) 
(14) 

(15) 
(16) 

The operator O5 ^ describes effects where one chirality flip occurs on the light quark side. Its matrix element scales 
like Qm s . 

There are no contributions scaling like Qm c , since the dependence on the charm quark mass must contain only 
even powers of m c . The leading contributions containing m c scale like m 2 and come from operators similar to @ and 
(fTTTfl . We will define them as 

= m 2 c s L [^-q^/q 2 ]h vL (17) 





-1) 


= rn b s L [iD ^- v^v ■ iD)]h vR 


ol- 


-1) 


= mb{v-id)s L [^^-v fl f/\h vL 


oi~ 


-1) 


= mbiid^ - v^v ■ id)]{s L h vR ) 


oi~ 


-1) 


= m b id u (s L [l^- v^i/\^ u h vB ) 


oi~ 


-1) 


= rn & m s s_R(7 M - v^i/)h vR , . 
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Operator Power counting Order in matching 

o[; 2 2) 1 q^(q) 

OtV A/Q a s (Q) 

Ofl ml/Q 2 a° B (Q) 

0% A 2 /Q 2 a°(Q) 

Of mj/Q i q°(g) 



TABLE I: Contributions to the long-distance amplitude for b — > s£ + £ coming from the different operators in the OPE Eq. |0J , 
together with the order in a s (Q) at which they appear in matching. 



C>2 0) = i m b^-s L <Ja V q v h vR . (18) 
<T 

There are many operators whose matrix elements scale like A 2 ; generally, they are of the form = 
qT(iD^)(iD u )h v or contain one factor of the gluon tensor field strength qTgG^ u h v . The latter operators can appear 
at O{a s ) in matching from graphs with qq quark loops as shown in Fig. 2(c), and can contribute to the B — ► K*£ + £~ 
amplitude through the graph in Fig. 1(a). 

Another class of operators appearing in the OPE describes effects of propagating charm quarks (see Fig. 1(b)), and 
have the form 

{2) = -^ I {sTh v ){cT c iD tl c). (19) 

The explicit form of these operators will be given in the next section, where it is shown that their contributions are 
further suppressed by m^/Q 4 relative to the short-distance amplitude. 

To sum up the discussion of this section, we argued that the long-distance effects to b — > s£ + £~ decays in the zero 
recoil region come from well-separated scales satisfying the hierarchy m& ~ Q > m c > A. These effects can be resolved 
using an OPE as shown in Eq. |JSJ. The contributions of the various operators in the OPE, relative to the dominant 
short-distance amplitude, are summarized in Table [H together with the order in matching (in a s (Q)) at which they 
start contributing. 

Some of the subleading operators appearing in the OPE give spectator type contributions to the exclusive B — > 
K*£ + £~ amplitude, as shown in Fig. 1. For example, the 0{AQ) operators o\ ^ and 0(A 2 ) operators can 
contribute through the graphs in Fig. 1(a), and the charm operators of the type Eq. lfH?fl contribute as in Fig. 1(b). 
Such spectator type contributions were studied at lowest order in perturbation theory in 0] where they were shown 
to be suppressed at least by A/Q. The effective theory approach used here extends this proof to all orders in a s , 
and shows that the suppression factor is a s {Q)A/Q (for the contributions from O)- 1 ') and A 2 /Q 2 (for contributions 
coming from O^ "*). 

We comment briefly on an alternative approach used in Refs. 0, ^| where the charm quarks and the large scales 
y/q 2 , nib are integrated out simultaneously. Such an approach includes the charm mass effects to all orders in m 2 /m 2 , 
but has the disadvantage of introducing potentially large power corrections ~ A 2 /m 2 . For this reason we prefer to 
integrate out only the large scale Q and leave the charm as a dynamical field in the OPE. 

The main result of our paper is that the contributions of leading order O(l) and the power suppressed terms 
0(m 2 /Q 2 ) to the long-distance amplitude depend only on known form factors and thus can be included without 
introducing any new hadronic uncertainty. The power suppressed terms of 0(A/Q) can be accounted for in terms of 
the form factors of the two dimension-4 currents qiD^^^hy. 

In the next Section we compute the matching conditions for these operators at lowest order in perturbation theory. 



III. MATCHING 

Typical lowest order diagrams contributing to the T— products T^(q 2 ) in QCD are shown in Fig. 1. The matching 
conditions for the operators appearing in the OPE Eq. JSJ are found by computing these graphs and expanding them 

in powers in 1/Q. At lowest order in a s (Q) the graph in Fig. 1(a) will match onto Oj , but not onto the O(AQ) 

operators Oj . These operators appear first at 0(a s {Q)) from graphs containing one additional gluon as shown in 
Fig. 1(b). 



6 








1 

2^ 




- tfiPLb) 


C#V)> = 


1 




- (ffiPLb) 
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FIG. 2: Graphs in QCD contributing to the matching onto sYh v operators (a), sYiD,j,h v (b) and sgG^T^hv operators (c). 
The filled circle denotes the insertion of Qi-e- In (c) the wavy line is the virtual photon 7* and the curly line denotes a gluon. 

An explicit computation of the graph in Fig. 1(a) with one insertion of the operators Qis gives the following results 
for the matrix elements of the T— products T i fJ '(q 2 ) on free quark states |8lll0|| (we use everywhere naive dimensional 
regularization (NDR) with an anticommuting 75 matrix) 

-| + | + 4G(m c )| (20) 

2 2 4 1 
-9^ + 9 + 3 G H (21) 

-^- e -l + ^G{m c )- 2 -G{Q)- % -G{m b )} (22) 

|_|_ 2G (0) + ! G (m c )-~G(m i) )} (23) 

_1 + AG{m c ) - 2G(m b )\ - -L m6 («(i% - q^P B b) (24) 

™ + \G{m c ) - r 3 G(m b )\ - ^m b (s(^ - q^)P R b) . (25) 



We denoted here with G(m q ) the function appearing in the basic fermion loop with mass m q l 



G{m q ) 



dxx{l — x) log 



-g 2 :r(l — x) + m 2 — it 



(26) 



In the kinematical region considered here (4m 2 < q 2 < 4m 2 ), this function is given explicitly by 



1 (m%\ 5 2m 2 



G(0) 
G(m b ) 

where in G(m c ) we denoted r = yl — Am 2 /q 2 



5 

18 



1, /XAJL 

6 ° g ^ 2 y 18 3q 2 



2m 2 



2to? 



1 1 + ^ • 
log — ^-«r 



1 /4m 2 



1 1 



2mf 



arctan 



4m~ 



(27) 
(28) 

(29) 



1 This function is related to h(z, s) used in ll(l as h(m q /m b , q 2 /m 2 ) = — 8/3G(m 9 ) — 4/9. 
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FIG. 3: Graphs contributing to the matching onto operators with explicit charm fields (see, e.g. Eq. 1331 1. In (a), (b) the 
filled circle denotes one of the QCD operators Qi_6- The crossed circle in (c) denotes the local operator appearing in the OPE 
with quark content (s&)(ce). The wavy line is the virtual photon 7* connecting to the e + e~ lepton pair. 



To match onto the operators introduced in Sec.|n]we expand the results (|20[1 - (|25[1 in 1/Q and go over to the HQET 
for the heavy quark field. To the order we work, this amounts to expanding the charm quark loop using 



G(m c ) = G(0) 



2 r 



log 



(30) 



On the other hand, since we treat ml and q 2 as being comparable, the full result for the b quark loop function G(m{,) 
has to be kept. 

To illustrate the matching computation we show how the result 1121 III for the T-product containing Qi is reproduced 
in the operator product expansion (jSJ). Expanding (|20() in powers of ml/q 2 one finds 



m = 



+4- 



log 



4G(0) 
1 



(31) 



O 



The terms of 0(q 2 ) and 0(ml) in this result can be identified with the matrix elements of the operators o[ ^ and 



respectively, provided that their Wilson coefficients are taken to be 



cfi a) 00 



1 

2^2 



4G(0) 



(32) 



Reproducing the 0{m^,/q A ) term in (|31|l requires the introduction of dimension-6 operators containing explicit factors 
of the charm quark field. They are obtained by matching from diagrams where the photon attaches to one of the 
external quark legs (see Fig. 3). Expanding these graphs in 1/Q and keeping only the term of 0(m c /Q 2 ) gives (the 
leading term scales like ~ 1/Q, but its b — > s matrix element vanishes) 



Tf(q 2 ) 



q- L 



(sjvP L b) 



(33) 



We dropped here operators which vanish by the equation of motion of the charm quark field (ip — m c )c = 0. The 
matrix element of this operator is computed by closing the charm loop, which gives 



(s\0 



r2^ b) = NcQc 



(<?V - q^P L b)-f { - + - - log ( -f 



1 



(34) 



The coefficient of the logarithmic term log m c agrees with that in the expansion of the exact result in Eq. I|31f) . This 
shows that the four-quark operators Eq. (|^3l reproduce the IR of the full theory result. However, these contributions 

are suppressed by m^/Q 4 ~ 0.8% relative to those of the leading operators O^', so they can be expected to be 
numerically small. This is fortunate, since their matrix elements on hadronic states would introduce new unknown 
form factors in addition to those contributing to the short-distance amplitude. In the following we will not include 
4-quark operators similar to those in Eq. (|33|l . 
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Using a similar expansion one finds the matching for all remaining T-products in H20(l - 125() onto the operators in 
the OPE ®. The results for the Wilson coefficients C\~ 2fl) are 



-n — C 2~i\v) 



C^i 2) (M) 



1 

2^2 

1 

1 

2^2 

1 

2^2 

1 

2^2 

1 

2^2 



4G(0) + - 

4 2 
>> + § 



(4tt) : 



f 36G( 9 2 ) . 
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G(0) 



1 

27 
7 



(4tt) 3 
- f G(m&) 



(-245(92) + 12C(g 2 )) , 



-^G(0)---|G(m„) 



4G(0) - 2G(m b ) - 



~G(0) - ^G(m b ) + ~ 



(35) 
(36) 

(37) 
(38) 
(39) 
(40) 



and 



CgV) = -^i 1 '^ 1 '^ 1 ^} (i = l-6) 
7r z 3 3 3 



(41) 



To facilitate the inclusion of the next-to- leading corrections, these results were computed using the operator basis in 
Ref. [23, and transformed to the basis in Eq. @ using 4-dimensional Fierz identities. For this reason, the constant 
terms in these expressions differ from those in Eqs. I|20(l . With this convention, the Wilson coefficients Gi(/z) used 
in the remainder of this paper differ beyond the LL approximation from those in Refs. [9J, Il0| and are equal to the 
"barred" coefficients G(/i) defined in Eq. (79) of Ref. (2i|. We included here also the next-to- leading results for 
C[~^ and C { 2 ~ 2 \ which can be extracted from the recent two-loop computation of SeidelyJ] (extending previous 



approximate results in 25]). The functions A(s), B(s), G(s) are given in Eqs. (29)-(31) of |24| and can be written as 



104 



243 



A(q 2 ) = -^log^ + 6A(q 2 ) 
B(q 2 ) 



Am 2 



243 ^ q 2 



243 



34 - 17™) log -4 



9 m t 

!1 °S 



1-7 1 « 2 , ml- 
17 log — log — f- 



C(q 2 ) 




SB(q 2 



16. q 
log — 

81 8 n 



The terms 8 A and <5-B do not contain explicit /i dependence and take the following values at the zero recoil point in 
B -> J ft:*e+e- (for (i = 4.8 GeV and m b (m b ) = 4.32 GeV) <5A(^ ax ) = 0.736 + 0.836i, (5B(g^ ax ) = -1.332 + 3.058«. 

Cq V \[i) is one of the Wilson coefficients appearing in the matching of the vector current (pj^b onto HQET currents 
and is defined in Eq. I|46|) . It accounts for the factorizable two-loop corrections not included in Ref. |24|. 
The results for the coefficients C\ 2 (n) can be computed in a similar way with the results 



(M) 



G<- 2) ( M ) = 0(a 2 ), 



Qia (m) 



(4tt) : 



-(48A(g 2 )), 



C 



(-2) 

i.a 



, , 12 2 11 

M = ^ { -9'"3'18'6 } (l 



(42) 



The Wilson coefficient (/x) appears in the matching of the tensor current qia^b onto HQET operators and is 
defined in Eq. (|4*T|) . The 0{a s {m b ) terms in the first two coefficients have been extracted from Ref. [2J|, where they 
are given in terms of the function A(q 2 ). 
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The only dimcnsion-4 operators appearing at this order in matching are 0\~ll and are introduced through the 
matching of the b field onto HQET according to b = (1 + -^~)h v . Their Wilson coefficients are 

CfcV*) = -C|- 2) (m b ) , C\- I] {m b ) = -C\f\m h ) = -C^ViO (i = 1 - 6) . (43) 

At two-loop order in the matching, all the other dimension-4 operators will appear, through the dependence of graphs 
such as those in Fig. 2(b) on external quark momenta. 

The gluonic penguin Q 8 contributes to the long-distance amplitude at leading order in 1/Q through one-loop graphs. 
The corresponding one-loop graphs were computed in the second reference of |25j in an expansion in q 2 /m 2 , and in 
Ref. [29j for arbitrary q 2 . Its contributions to the Wilson coefficients of the leading operators are 

(tx 2 V) - jf^i'V) , cifM = -ggltf) (q 2 ) (44) 

with F 8 (9 ' 7) (g 2 ) given in Eqs. (82), (83) of Ref. |2^. The operator Q 8 contributes also at tree level through gluon- 
photon scattering graphs (with the photon coupling to the b and s quarks). Expanding these graphs in powers of 1/Q 
one finds at leading order 

M ^ - mbQfc s L a af3 gG a ^h vR + ^SLl^a a pgG^h vR . (45) 
(47r) z w • q 8ir^ 

The matrix elements of these dimension- 6 operators are suppressed by A 2 /Q 2 . 

The one-loop graphs in Fig. 2(c) with one insertion of Qis produce dimension-5 operators containing the gluon 
field tensor of the form sgG^h v . Although their Wilson coefficients start at 0(a°), their matrix elements are ~ A 2 , 
and therefore are suppressed by A 2 /Q 2 relative to the short distance amplitude. We will neglect all these higher 
dimensional operators and keep only the O(l), 0(m 2 /m 2 ) and 0(A/rrib) terms in the long distance amplitude. 



IV. MATRIX ELEMENTS 



In this section we use the OPE result Eq. (jHJ for the long-distance amplitudes Ti(q 2 ) to compute the hadronic 
amplitude in Eq. J3J up to and including corrections of order 0(a s (Q) 1 K/mb 1 rn 2 /m 2 ). At this point we 

encounter a technical complication connected with the fact that the OPE was performed in terms of HQET operators, 
while the matrix elements of the QCD currents sTb appearing in the factorizable matrix elements of Qr t g are expressed 

(—2) 

in terms of physical form factors. This means that the matrix elements of the operators 0{ 2 are given in terms 
of HQET form factors, which are not known. Also, keeping all 0(A/m&) contributions requires that we include also 

( — 2) 

T-products of the 0\ 2 operators with l/rn& sub-leading terms in the HQET Lagrangian. Such nonlocal matrix 
elements introduce additional unknown form factors. This proliferation of unknown matrix elements appears to 
preclude a simple form for our final result. 

We will show next that it is possible to absorb all these nonlocal matrix elements into the physical form factors, 
through a simple reorganization of the operator expansion, such that one is left only with local 1/mb corrections. 

This can be achieved by expressing the leading operators 0\ in terms of QCD operators, up to dimension-4 HQET 
operators siD^ipf^hy. Technically, this is obtained by inverting the HQET matching relations (we assume here 
everywhere the NDR scheme) 

slI^l = C^ v \fx)s L ^^h V L + C^\^)s L v^h vR + -— SLJuiPhvR + 0(1 /ml) (46) 

SLicr^q v b R = CQ t \fx)sLi(T^uq v h V R + C[ t \fi)SL[(v ■ q)^^- ^Vfj,]h V L (47) 

+ -^—s L ia l _ lv q u iph vL + 0(1/ 'm 2 ) . 
2nib 

The Wilson coefficients C^' 1 ' (p) are given at one-loop by [2(J 

C#V) = 1 
Cf (/x) = 1 



— 3 log — +4 , 

47r nib 

— (5 log — +4), 

4-7T mb 



Ust>F_ 

2tt 



Cf = 0(a 2 s ) . 



(48) 
(49) 
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In the 0(l/rrib) terms we work at tree level in the matching, which will be sufficient for the precision required here, 
although the method can be extended to any order in a s (mj,). 

Solving the matching relations Eqs. I)4(i|l. 147(1 for the leading order HQET operators appearing in the OPE o[ ^ 
one finds 



sl(9 2 7/x - <i<ln)hvL 



SLia^q^hyR 



1 



1 C^\fi) 



1 



(-1) 



1 



-o 



(-1) 



(50) 
(51) 



We neglected here terms of 0(a s (mb)A/mb)- 

Substituting these results into the OPE, the leading terms can be written in terms of physical B — > K* form factors, 
with corrections of 0(A/rrib) coming from local dimension-4 operators o[_^ 



A (V) 



-Cf(t,)^(sia^(l+ l5 )b) 

1 5 



(52) 



C-i)\ 



We absorbed here the contributions from the leading terms in Eqs. I]50p. (|51|l into a redefinition of the Wilson 
coefficients C7 9 



6.8 



c 7 (m) 
c 9 (aO 



Cf(/x)=C 7 ( )U ) + 27r 2 ^C i ( M ) 

i=l 

Cf( M ) = C 9 (jm) 

6,8 



(53) 
(54) 



The O(l) and 0(m 2 /q 2 ) contributions to the long-distance amplitude are contained in C% 9 , and the 0(A/m b ) part is 
encoded in the matrix elements of C)| . Note that the effective Wilson coefficients Gy* introduced here are different 
from the "effective Wilson coefficients" commonly used in the literature Cj% 0, [25| . The latter include contributions 
from the matrix elements of the operators <5i-g (usually computed in perturbation theory), and are thus dependent 
on the final state. In contrast, our effective Wilson coefficients are state independent, and encode only contributions 
from the hard scale /1 ~ mj. 

Combining everything, the next-to-leading expressions for the effective Wilson coefficients are 



Cf = C7 9 -(C 1 + ^)[8G(0) + |]-C 8 [f G(0)-yGH + A] + C 4 [^(0) + yGK) + ^] 



(55) 



G 5 [8G(0) - 4G(m fc ) - ^] - G 6 [^G(0) - ^G(m b ) + ^} + ^-[ Cl 9C(q 2 ) + C 2 (-6B(q 2 ) + 3C(q 2 )) - G 8 ifV)] 



off 



G 7 - ^G 3 - ^G 4 + l -C 5 + l -C 6 + ^-[-C 2 6A(q 2 ) - G 8 if V) 



(56) 



The effective Wilson coefficient G| ff is RG invariant. At the order we work here, it satisfies the RG equation 



f i-^Cf(n)=0(a 2 s C 1 , 2> a s C 3 ^). 



The coefficient G£ ff (/i) satisfies a RG equation 

d 



M^G 7 cff ( Ai )= 77 (a s )Gf( M ) 



(57) 



(58) 



with anomalous dimension 77(a s ) = 7*(a s ) — 7 m (a s ) (see Eqs. (|A32() and (|A33(1 for definitions). 
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The Wilson coefficients of the dimension-4 operators Bi (p) are given by 

BiO*) = 8 7 r 2 ^Q( M )(C| > 7 1) (M)+c| i 7 2) ( M )) (59) 

i 

B 2 (p) = Stt^QMC^V) (60) 

i 

B 3 Qj) = S^a^C^V) (61) 

i 

B^) = s^^Ci^J^^-icr^V)) (62) 
S 5 („) = STr^^^^^ + i^^)). (63) 

i 

These Wilson coefficients start at 0(a s ) in matching. By absorbing the factor of Sir 2 in their definition, their expansion 
in a s (Q) starts with a term of order a s (Q)/TT. At the order we work (keeping terms in the OPE of 0(a s , A/Q, m 2 /Q 2 ), 
but neglecting 0(a s A/Q) terms), they all vanish Bi_5 = 0. However, we will include them in the following expressions, 
which is required for a complete result to 0(A/rrib) accuracy for the long-distance amplitude. 

It is convenient to parameterize the physical amplitudes A^'^ introduced in Eq. in terms of eight scalar form 
factors A {V > A) - V (V ' A ) defined as 

A (V.A) = A (VA) iq 2 )t£ ^ XaV *, {p + k) X {p _ kr + B (V,A) {q 2 ) ^ ( g 4) 
+C^ A \q 2 M ■ p)(p + k), + V^ A \q 2 M -p){p~ h)» 

The B — > K*e + e~ decay rate can be represented as a sum over the helicity A = ±1, of the vector meson. In the 
limit of massless leptons, this is given by 



dT(B -> K*e+e-) _ 4G 2 F \V tb V t *\ 2 a 2 



q 2 \q\ £ {\H { ^\ 2 + \Hi A) \ 2 } (65) 



dq 2 3to 2 b (4tt) 5 

H BK ' A=±1,0 

where the i?\ and correspond to the vector and axial leptons coupling, respectively. Expressed in terms of 

the scalar amplitudes A, B,C introduced in Eq. I)64[l. they are given by (i = V, A) 

H { l\q 2 ) = T2m B \q\A^{q 2 )-B^{q 2 ) (66) 
i^V) = -^^{{-q 2 +ml-m 2 v )B^{q 2 )+A m 2 B q 2 C^{q 2 )X . (67) 



2mv\/q 

The explicit results for the amplitudes A^ V,A ^ — V^ V ' A ^ are obtained by taking matrix elements on physical states 
and are given by 



A^(q 2 ) 


q A 


+ (q 2 


B^ V \q 2 ) 


cf(») 2m 2 b i( 




C^ V \q 2 ) 


q A 


-<?+( 


V^ V \q 2 ) 


= 0*)— a-Iff- 


(q 2 ) 



-reff / 



-.off/ 



-rcff / 



(68) 
(69) 
(70) 
(71) 



and 



A (A Hq 2 ) = C 10 g(q 2 ) (72) 

B^(q 2 ) = -C w f{q 2 ) (73) 

C {A) {q 2 ) = ~C w a + {q 2 ) (74) 

V^(q 2 ) = -C w a_{q 2 ). (75) 
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The coefficients V^ V,A ^ (q 2 ) do not contribute to the B — > Ve + e decay rate into massless leptons, but are relevant 
for the B — > Vt + t~ mode. We will not consider them further. The 0(A/Q) contribution to the long-distance 

contribution appears as matrix elements of the local dimension-4 operators 0\ (denoted as A\.d.{q 2 ) — 2?i.d.(<7 2 ) in 
Eqs. (jHEJ-lZIJ)- They are given explicitly by 

Ad.Oz 2 ) = ^-d^{^)B 1 + ^-{A-v.k)g{q 2 )B 2 -^-[{l + A.) g+ tf) + (i- A) ff _ (g 2 )]i?4 (76) 

ZTOft iTOf, 4771;, m B m B 



2mb 



9(q 2 )B 5 



Bi. d .(q 2 ) = 7 ^-df , \q 2 )B 1 --^-(A~v.k)f(q 2 )B 2 (77) 
2 nib ^ m b 

1 _ rffi 

[(A« ~ k).(p + k)g + (q 2 ) + (Av - k).(p - k) 9 -(q 2 )}B 4 + f^f(q 2 )B 5 

Cl.d.(g 2 ) = ^[dfiq^-^^siq^-^A-v.k^a^ + ^-siq^ (78) 
Zrrib im,b Irrib imb 



1 , v.k „ , t n „ , A , ON 1 /T , \ 



-(1 - — )s(q 2 )B 3 + (-^ 9+ (q 2 ) - -(A - v.k)h(q 2 ) - — s{q 2 ))B A 
4rn{, m B Ami I Ami, 



a+ (^) + — s(^)]B 5 



2mfc 2wb 

The form factors appearing here are defined in the Appendix. The corresponding result for £>i.d.(<7 2 ) can be obtained 
from the Ward identity which gives £>i.d.(<7 2 ) + ( m s — m v)C\.d.{q 2 ) + q 2 T^\.d.(q 2 ) = 0- Expanding in powers of 1/rrib 
and keeping the leading terms gives Z\d.(<7 2 ) = — Ci.d.(? 2 ) + 0(A/mb). 

For completeness we quote here also the relevant results for the semileptonic decay B — > pev. The decay rate is 
given by a sum over contributions corresponding to helicities of the final vector meson A = ±, 



dT(B^peis) _ G 2 F \V ub \ 2 2 



q q 



E l^l 2 (79) 



dq 2 967t 3 TOr 
* B A=±1,0 

where the helicity amplitudes are given by 

H ± {q 2 ) = T 2m B \q]g( q 2 ) + f(q 2 ) (80) 

H (q 2 ) = - 1 {(g 2 - m| + m 2 v )f(q 2 ) - Am 2 B q 2 a+{q 2 )} . (81) 
2myy q 



PHENOMENOLOGY 



In the low recoil region, the amplitudes A^ v \ B^ v \ for B — > K*£ + £~ are dominated by the operator Qg. 
The contribution proportional to CV can be expressed in terms of the Cg terms using the form factor relations 
Eqs. ljA23() . (|A38|) . \ A39(l given in the Appendix. Keeping terms to subleading order in A/ nib, these amplitudes can be 
written as 

AWtf) = Cfg(q 2 ) 
B^(q 2 ) = ~Cff(q 2 ) 
C^(q 2 ) = 



Cfa+(q 2 ) 




Inserting these results into the expressions for the helicity amplitudes H^\q 2 ) one finds 
H { P{q 2 ) = T2m B m v ^F^lCfg{q 2 ){l + 5 + r a ) + Cff{q 2 ){l + 5 + rb) 



<V) = - 



m B y 



-Cff(q 2 ){l + S + r b )- 2m 2 B m v - 



■Cfa + (q 2 )(l + 6 + r c ) 



(85) 
(86) 
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Here 1 + S(q 2 ) scales like m° and r abtC (q 2 ) parameterize the l/m b correction. Their explicit expressions are 



-<ef 
y 9 



^(g 2 ) 


= % 2 )— ( 


r b (<z 2 ) 


^ % 2 ) — ( 

m b 




= % 2 ) — ( 
m b 




(88) 
(89) 

Cf a+(g») ■ (90) 

Combining the RG equations satisfied by Cf ff (/u) Eq. and by the 1 + 2D^\p)/Cq V \/j) factor Eq. <|A31|) . one can 
see that the S(q 2 ) parameter is RG invariant. 

These results imply that the H^\q 2 ) amplitudes for rare B — > V£ + £~ decays are related at leading order in A/m& 
to those for semileptonic decay B — ► Vev with a common proportionality factor 

#f V) = Cf (1 + 5{q 2 ) + 0(A/m b ))H X (q 2 ) . (91) 

Combining this with the rate formulas one finds a relation among the decay rates for the rare and semileptonic decays 

dr(B-,pev)/dq 2 \V ub \ 2 St: 2 1 EaI<^V)! 2 (92) 



dT(B^K*e+l-)/dq 2 \V tb V t * s \ 2 a 2 \Cf{l + 5{q 2 ))\ 2 + \C w \ 2 Y,x\H^ K \q 2 )\ 2 

The corrections to this relation are of order 0(A/m b ) and can be expressed in terms of the three parameters r a btC (q 2 ) 
introduced above. 

The ratio of decay rates in Eq. I|92|l has been considered previously in Ref. [T3, Il3l Il4| in connection with a method 
for determining | V ub \. This requires some information about the SU(3) breaking ratio of helicity amplitudes appearing 
on the right-hand side 

Rb{v) = ET^^W (93) 

It has been proposed in |l2l|l4 to determine Rb in terms of the corresponding ratio of D — > p/K* decay amplitudes 
Rd{u) using a double ratio |15|. up to corrections linear in both heavy quark and SU(3) symmetry breaking 

R B (y) = R D (y)(l + 0(m s (- -)) (94) 

m c m b 

In this relation, the two sides must be taken at the same value of the kinematical variable y = Ey / my . A chiral 
perturbation theory computation |l2j | at the zero recoil point y — 1 shows that the corrections to this prediction are 
even smaller than suggested by the naive dimensional estimate Eq. I|94|l . We do not have anything new to add on this 
point, and focus instead on the structure of the denominator in Eq. I|92l) . 

The results of our paper improve on previous work in two main respects. First, we point out that the rate ratio 
(|92ll can be computed at leading order in l/m b over the entire small recoil region, and not only at the zero recoil 
point q 2 = (rriB — my) 2 - This has important experimental implications, as the rate itself vanishes at the zero recoil 
point, such that measuring the ratio in Eq. I|92|) would involve an extrapolation from q 2 < qf nliK . Most importantly, 
Eq. (|92|l allows the determination of V ub using ratios of rates integrated over a range in q 2 , as long as such a range is 
still contained within the low recoil region. 

Second, we present explicit results for the subleading 0(A/m b ) correction to this result in terms of new form factors 
contained in the parameters r atbiC (q 2 ). Using model computations of these form factors, this allows a quantitative 
estimate of the power corrections effect on the V ub determination. 

In the rest of this section we will study in some detail the (RG- invariant) quantity N c g(q 2 ) defined through the 
ratio of rare radiative and semileptonic decays in Eq. I|92() 

dr(B-,pev)/dq 2 \V ub \ 2 8n 2 1 £ A \H^{q 2 )\ 2 



dT(B^K*£+£-)/dq 2 \V tb V t * s \ 2 a 2 N eS (q 2 ) £ A \H* -** (q 2 )\ 2 
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fib (GeV) 


c 9 


c 7 


Cf(y = l) Cf(y = 1.5) 


Cf(y = l) 


Cf (y =1.5) 


N cS {y = 


1) N cB (y = 1.5) 




2.4 


4.378 


-0.388 


4.315 + 0.198i 4.338 + 0.198i 


cv 


C 7 


30.80 


28.96 


LL 


4.8 


4.140 


-0.343 


4.331 + 0.550* 4.395 + 0.550i 


cv 


C 7 


33.37 


32.34 




9.6 


3.760 


-0.304 


4.420 + 0.822* 4.513 + 0.822* 


cv 


c 7 


35.81 


35.38 




2.4 


4.510 


-0.366 


4.685 + 0.494* 4.742 + 0.442i 


-0.352-0.127i 


-0.360-0. 122i 


32.75 


30.83 


NLL 


4.8 


4.218 


-0.332 


4.611 + 0.556* 4.680 + 0.514i 


-0.401-0.100i 


-0.408-0.097i 


32.76 


31.11 




9.6 


3.799 


-0.300 


4.589 + 0.643* 4.668 + 0.609*' 


-0.422-0.083i 


-0.428-0.080i 


33.46 


32.10 



TABLE II: Results for the Wilson coefficients in the weak Hamiltonian CV,g and the effective Wilson coefficients appearing in 
the B -> K*e+e~ decay rate at LL and NLL order. The Wilson coefficient Cio is equal to Cf LL = -4.409 and C 1 N NLL = -4.279. 
The other parameters used here are m^ijni,) = 4.32 GeV, a a (Mz) = 0.119 and m c (m c ) = 1.335 GeV. 



The results of this paper offer a systematic way of computing this quantity in an expansion in a s (Q), m 2 /Q 2 and 
A/rUb- The precision of a \V u b\ determination using this method is ultimately determined by the precision in our 
knowledge of this parameter. There are several sources of uncertainty in N c ff(q 2 ), coming from scale dependence, 
0(A/rrib) power corrections and duality violations. We will consider them in turn. 
At leading order in A/m,b, the N e g(q 2 ) parameter is given by 



N eS (r) = \c e 9 + 



2m b {nf 



cs 



\ 2 + \C 10 \ 2 + O(A/m b ). 



(96) 



We give in Table II results for the effective Wilson coefficients C 7 9 at several values of the renormalization scale 
/i ~ m;, . We work both at leading log order (next-to-leading log order for Cg(/x)), and at next-to- leading order (NNLL 
order for Cg). In each of these approximations the combination of effective Wilson coefficients in Eq. (|96(l satisfies 
the RG equation 



d_ 

dfi. 



s 2m b (/i) 2 ^ eff 

°9 ' 9 °7 



2°° 



(v) 



(M) 



0(a s ), 
0(a 2 C h2: a s Q 



3-6) 



(LL) 
(NLL) 



(97) 



The structure of the NNLL running for the Wilson coefficients in the b — > se + e~ weak Hamiltonian was given in 
Ref. [23| (see also [2^ ). The complete NNLL result requires the 3-loop mixing of the four-quark operators into Q7,g, 
which was obtained only recently [3(j- We use here the full NNLL results for the Wilson coefficients , which were 

presented in [23|. The factor containing D^\fi) can be extracted from Eq. (|A24Jl and its inclusion is necessary at 
NLL to achieve the scale independence of N e g to this order. 

To illustrate the q 2 dependence of the effective Wilson coefficients, we quote their values at two kinematical points 
y = 1 and y = 1.5, corresponding to the low recoil region overlapping with that kinematically accessible in D decays. 
The resulting dependence on y is very mild, of about 2.5% in Cg S and almost negligible in Cj S . 

Next we consider the scale dependence of the results, by computing the variation of the effective Wilson coefficients 
between the scales 2/i;, and /i&/2 with /i^ = 4.8 GeV. The LLO Wilson coefficient Cg changes in this range by 15%, 
while the corresponding variation in Cg ff is reduced to 2% (for the real part), and 36% (for ilm Cg ). At NNLL 
the change in Cg is 17%, which is reduced in the effective Wilson coefficient Cg to 2% for Im C| , and 8.5% for 
ilm Cg . Combining everything, at LL order the scale dependence of N e g is about 16% which is reduced at NNL 
order to about 3.5% (at the zero recoil point y = 1). 

To get a sense for the relative contributions to the long-distance effects in C$ s , we give below the detailed structure 
of this effective coefficient at LL and NLL orders for fib = 4.8 GeV at y = 1 



LL : C 9 off (y 
NLL : C 9 off (y 



1) = 4.140 + (0.136 + 0.506i) + (0.004 + 0.044*) + 0.000 + 0.050 = 4.330 + 0.550* 
1) = 4.218 + (0.313 + 0.505*:) + (0.001 + 0.050*) - 0.006 + 0.085 = 4.611 + 0.556* . 



The five terms correspond to Cg, the contribution of CJ1.2, from Q3-6) Qs and the m 2 /Q 2 term respectively. As 
expected, the dominant contribution to the long-distance part of Cg ff comes from the operators <3i,2, with Qz-q 
contributing about 3% and the m 2 /Q 2 term about 0.1%. 

The structure of the power corrections of 0{A/mb) is in general very complicated and depends on both the leading 
and subleading B — > V form factors. Details of such an analysis will be presented elsewhere. We will limit ourselves 
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here to the study of these corrections at the zero recoil point, where they are given only by ri,(q 2 ), defined in Eq. I|89[l . 
At the zero recoil point q 2 = ^ lax , the relation among rare radiative and semileptonic hclicity amplitudes Eq. 191|) 
can be extended to subleading order in 1/ mj and reads 

h ( P (ql ax ) = cf (i + <%LJ + n.(<?Lx))ffA(<?LJ • (99) 

The corresponding modification of the relation for decay rates Eq. (|92[) is obtained by the replacement 1 + 5(q 2 ) — ► 
1 + S(q 2 ) +ri,(q 2 ). Since the leading order result for N e g(q 2 ) has only a weak dependence on q 2 in the low recoil region 
(see Table II), this is a reasonably good approximation. 

A complete computation of r b (q 2 iaK ) is not possible at present as Bj.d. depends on the (as yet unknown) Wilson 
coefficients Bi_$. Dimensional analysis estimates of the first term in l|89|l give ^(q^ax) ~ — (0.03± 0.01)A/m&, which 
represents at most an uncertainty of 1% in N c s (g^ ajt ). Barring an anomalously large value for B\.d., this suggests very 
small power corrections to the coefficient N c g. 

Finally, we address the issue of duality violations. Their effects are difficult to quantify in a precise way, but some 
guidance can be obtained from the experimental data on the R = a(e + e~ — > hadrons) / a (e + e~ — > u + u~) ratio, to 
which the coefficient N c g(q 2 ) is very similar. Good data is available on the ratio R in the cc resonance region (see 
e.g. Fig. 39.8 in j23). In the region y/q 2 = 4.1 — 4.4 GeV (corresponding to the kinematics relevant here), the ratio 
R oscillates around its pQCD predicted value by less than ~ 25%. Strictly speaking, the quantity analogous to R in 
our case is Im(Cg ff ), which represents only about 12% of the magnitude of \Cf |. In the real part of Cf, the relative 
error introduced by these oscillations is suppressed by the large value of Cg to about 0.3/4.3 x 10% ~ 1%. Due to the 
fact that Im (Cf)/Re (Cf) ~ 12%, the 25% duality violation effect in Im(C| ff ) is reduced in \Cf + 2m 2 /q 2 Cf\ 2 
to about 2%. The corresponding effect in N e g is reduced by a further factor of 0.5 since the contributions of the two 
terms in N e s are roughly equal, and Cio is an invariant. These arguments show that duality violation effects are 
likely to be very small in N e g in the kinematical region considered, probably below 5%. Precise measurements of the 
q 2 spectrum in this region could help resolve and reduce this source of uncertainty. 

Combining all sources of errors, we find a total uncertainty in N e g of less than ~ 10%, which is dominated by 
duality violation effects. This gives a total theory uncertainty on \ V u b\ from this method of about 5%. 

We comment briefly on the experimental feasibility of this method. Model estimates of the dilepton invariant 
mass spectrum in B — > K*/i + fi" indicate that the integrated branching ratio corresponding to the region considered 
here q 2 = [15, 19] GeV 2 is about (2 — 5) x 10~ 7 , depending on the form factor models used 0. Extrapolating the 
uncertainties in the present data 0,0] to 1000 fb _1 , corresponding to the entire data sample from the B factories, 
suggests that this integrated branching ratio will be measured to about 25%. This is beginning to be comparable to 
the theory uncertainty, and indicates that a competitive determination of \V u b\ using this method will likely require 
a super B-factory. 

VI. CONCLUSIONS 

We presented in this paper a short-distance expansion for the long-distance contributions to exclusive B — ► K^£ + £~ 
decays in the small recoil region. The main observation is that in this kinematical region, there are 3 relevant energy 
scales: Q — mj ~ \fq 2 , m c ,A. We use an operator product expansion (OPE) and the heavy quark effective theory 
(HQET) to integrate out the effects of the large scale Q, and classify the effects from the remaining scales in terms of 
operators contributing at a given order in m 2 /Q 2 and A/Q. 

Our main result is a systematic expansion for the long-distance amplitude in B — > K^*H + £~ decays including terms 
of 0(m 2 /Q 2 ) and 0(A/Q), which can be extended to any order in a s (Q). The final results for physical observables 
are explicitly scale and scheme independent, order by order in perturbation theory. This is to be contrasted with 
the often used naive factorization approximation (combined with resonance saturation), which is not consistent with 
constraints imposed by renormalization group evolution. 

The form of the result is analogous to that for the R ratio in e + e~ — > hadrons, which can be computed systematically 
in an expansion in 1/Q 2 . For example, the nonperturbative effects in the R ratio have an analog in the b —> se + e~ case 
as form factors of higher dimensional flavor-changing currents. We classify all the nonperturbative matrix elements 
required for a complete description of B — > K^*H + £~ to the order considered. We find that none of these new 
form factors enter at order O(l) and 0(m 2 /m 2 ) for the long-distance contribution, and start contributing first at 
0(a s (Q)A/m h ). 

These results are applied to a method for extracting the CKM matrix element V u b from the ratio of semileptonic 
and rare exclusive B decays in the small recoil region. We find that the long-distance effect in this determination 
is well controlled by the expansion in A/m& and m 2 /m 2 , and the precision of such a method is dominated by scale 
dependence and duality violating effects. Experimental measurements of the dilepton invariant mass spectrum dT/dq 2 
in B — > K^ > f r t~ will allow a direct control of these effects. 
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The methods of this paper can be applied to other problems of interest for the phenomenology of rare B decays. 
The long distance amplitude has a complex phase, which is however completely calculable using the OPE. This means 
that observables such as CP violating asymmetries (in the SM and beyond) can be computed in a model-independent 
way. Combined with methods based on the soft-collinear effective theory (SCET) |3^| and perturbative QCD |29l l34| , 
which are applicable at large recoil, the approach proposed here opens up the possibility of attacking the exclusive 
b — > se + e~ rare B decays from the both ends of the q 2 spectrum. 
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APPENDIX A: FORM FACTOR RELATIONS 



We give here an alternative derivation of the improved heavy quark symmetry form factor relations at low recoil 
presented in Ref. |||, including the leading power corrections ~ 0(A/m,£,) and hard gluon effects. As a by-product 
we derive exact relations for the HQET Wilson coefficients of dimension-4 operators following from the equations of 
motion. 

We start by giving the definitions of the B — > V form factors used. One possible parameterization is 



(V(k, v )\q^b\B(p)) = g^ie^x^ip + kfip-ky 

ji \ /„* 



(V(k,r,)\qw 5 b\B(p)) = f [q 2 )^ + a+tfW ■ p){p + k) ^ 

+ a_{q 2 )( V * -p)(p-fc) M , 
(V(k, r))\qia^b\B(p)) = g + (q 2 )ie^ XaV * x (p + kf + g_ (q 2 )ie^ Xa rf x (p - kf 
+ h(q 2 )(rf ■ p)ie^xa(p + k) x (p - k)° . 



(Al) 
(A2) 

(A3) 



We use the convention £ 0123 — 1. This particular definition of the form factors is convenient in the low recoil 
region q 2 ~ (tob — mv) 2 , where it simplifies the power counting in mf,. Taking into account the usual relativistic 
normalization of the B meson state, these form factors satisfy the scaling laws Q 

tl 2\ 1/2 / 2\ -1/2 , 2\ ( 2\ -1/2 i 2\ , ( 2\ -3/2 

J{q ) oc m b ' , g{q ) oc m b a+{q ) - a^(q ) oc m b , a+(q ) + a_ {q ) oc m b 

g+{q 2 ) - g-(q 2 ) «■ m l /2 > g+{q 2 ) + g-(q 2 ) « m b 1/2 > h (q 2 ) <* m b 3/2 ■ ( A4 ) 

We will require also the form factor of the pseudoscalar density defined as 

(V(k, V )\q l5 b\B(p)) = ( V *-pHq 2 ). (A5) 

This is not independent and can be obtained using the equation of motion for the quark fields in terms of the form 
factors defined above as 



S(q 2 



1 



m h + m q 



{f(q 2 ) + (m 2 B - m 2 v )a+(q 2 ) + q 2 a„{q 2 )} . 



The leading term in the expansion of s(q 2 ) in powers of A/rrib scales like s(q 2 ) oc m b 1 and can be written as 

■<q 2 ) = -f(q 2 ) - a+(q 2 ){m B + v ■ k) - a-(q 2 )(m B - v ■ k) + 0(m b 3/2 ) . 

m,B 



(A6) 



(A7) 



An alternative parameterization commonly used in the literature defines the form factors as (with gy, = p^ — fc„) 

(A8) 



(V(k,r,)\q^(l-~/ 5 )b\B(p)) = 2V(q2) ie^ P *V w P p k° 



TUB + my 

-2m v A (q 2 ) q» - {m B + m v ) A 1 (q 2 ) 



+ A 2 (q ) ■ 



„1 Qf 1 



TUB + my 

{Vik^qia^q^Bip)) = -2T 1 (g 2 ) ie^ p(y r, w p? k° 
{V(k, vM^uq^blBip)) = T 2 (q 2 )[(m 2 - tt^Jtj* - (rj* • p)(Pp + *:„)] 

+ T ^{q 2 )—T — ^r[( m B - m v){p^ - M - q 2 (Pn + M] 

m B m V 

The relation to the alternative definition in Eqs. l)Aljl - l)A3|l is 



(A9) 
(A10) 



g(q 2 ) = - 

a + {q 2 ) - 



1 



tub + my 



1 



m B + m v 
g + uf) - T^q 2 ), g (q 



V(q 2 ), f(q 2 ) = (m B +m v )A 1 (q 2 ) 

2 2 2m V 2 TOi3 + my 2 m B - my 
^ 2 (g ) , a-{q ) = -^-A (q ) 5 A nQ ) H 2 



(All) 



^2 (<Z 2 ) 



2 \ rp /2\ „ /„2\ _ m B m V /rp (2\ rp I „2\\ L/„2\ _ ^ //Tl /„2\ /„2-l 



-(T 2 fy) - r^q 2 )) , = -3iTi{q 2 ) - T 2 (^)) 



m B — m 



-Uq 2 ). 



B _ "V 
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In addition to these form factors, we require also the matrix elements of the dimension-4 operators qiD ^(75)6, 
which can be defined as 

(V(k, r])\qiD fl b\B(v)) = d(q 2 )ie ^rf* '(p + fc) A (p - kf (A12) 
(V(k,r ) )\qi < D fi j 5 b\B(v)) = d^q 2 )^ + d+(q 2 )( V * ■ p)(p M + fc„) (A13) 

+ d-(q 2 )(jf -p)(p p - fc M ) . 

Their scaling with the heavy quark mass is complicated by the presence of the covariant derivative i-D^, which can 
introduce factors of the large scale m b through loops. To make it explicit, we consider the matching of the dimension-4 
QCD operators in Eqs. I|A12|) . (|A13(1 onto HQET. Working at tree level in the dimension-4 operators, but keeping all 
contributions enhanced by 0(rtib), this can be written as 

qiD^b = D ( Q v \p)m b q^^h v + D[ v \fi)m b qv^h v + qiD ^h v H (A14) 

qiD^ab = -Do (AO m i> 97^75^ + D^\p)m b qv^ 5 h v + qiD ^ b h v H (A15) 

We assumed here the naive anticommuting 75 scheme. The Wilson coefficients d[ v \[i) start at 0{a s ). 

The matrix elements of the dimension-4 HQET operators analogous to those appearing in Eqs. (|A12|I . (|A13(I 
(obtained by replacing qiD^(j 5 )b — ► qiD ^(7s)/i„) can be parameterized in terms of similar form factors, denoted 
with d,(°\q 2 ), .... They have a simple scaling with the heavy quark mass, which is the same as in Eq. (|A4|) with the 

substitution (d( ),4 0) ,4 0) > d - 0) ) -> ($»/> 

a+, a_). These form factors are related to the effective theory form factors 

introduced in [6( as 

d^(q 2 )= 1 -V(q 2 ), 4 V) = -2M<? 2 ), ••• (A16) 
Taking the B — ► V matrix element of Eq. (|A14|) one finds for the leading terms in the 1/rrib expansion of d(q 2 ) 

d(q 2 ) = ^^-m b9 (q 2 ) + d<°> (q 2 , M ) + . . . (A17) 
Co (W 

We keep here all terms of order 0(a s m^ 2 ) and 0(m b 1 ^ 2 ) and the ellipses denote terms of order 0(a s m b 1 ^ 2 , m b 3 ^ 2 ). 
Similar expansions are obtained from Eq. (|A15(1 

d ^ 2 ) = - £ j^rn b f(q 2 ) + df\q 2 ) + ... (A18) 

= -^\^ + {q 2 ) + df{q 2 ) + ... (A19) 
Co (M) 

= -^^^(q 2 ) + d^(q 2 ) + ... . (A20) 

In the low recoil region, heavy quark symmetry predicts relations among these form factors 0, ■ The sub-leading 
corrections to these relations were computed in We give here an alternative simpler derivation, valid to all orders 
in l/m b (see also [I^). We take this opportunity to include also 0(m q ) light quark mass effects (with m q the mass 
of the quark produced in the weak decay b — > q) in these relations, which were neglected in Such effects can be 
important for the case of B — > K* decays. 

The first relation is obtained from the operator identity 

id v {qia liv b) = -(m b + m q )qj,j,b - 2qiD ^b + id^qb) , (A21) 

which follows from a simple application of the QCD equations of motion for the quark fields. Taking the B — > V 
matrix element one finds the exact relation 

g+(q 2 ) = -{m b + m q )g{q 2 ) - 2d{q 2 ) . (A22) 
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Counting powers of nib and keeping the leading order terms gives the well-known Isgur-Wise relation among vector 
and tensor form factors Q g+(q 2 ) = —m,Bg{q 2 )- Keeping also the subleading terms of 0(m b 1 ^ 2 ) reproduces the 
improved form factor relations derived in Inserting the expansion of d(q 2 ) Eq. (| Al 7|> into Eq. (|A22Jl gives 

9+(q 2 ) = ~ U + S^^l "W) - M( % 2 ) ~ m ^(q 2 ) + ■■■■ (A23) 
V Co {n) J 

This agrees with the improved symmetry relation Eq. (48) of Ref. Q and generalizes it by including light quark mass 
effects and by making explicit the renormalization scale dependence. The radiative corrections to this relation were 
computed in Ref. @] at fi = mb in terms of a coefficient k,\ (defined in Eq. (23) of |(|). Using Eq. (|A25|) below this 
coefficient can be expressed as 

The equation of motion Eq. 1A21|) can be used to determine the Wilson coefficients Dq\ (/i) in the matching of the 
dimension-4 operators Eq. (|A14|) in terms of the Wilson coefficients of the dimension-3 currents. In this derivation 
we set id^ = nisv^ — = niBV^(\ + 0(A/m{,)). We find 



Cf (M ) - Cf 0*) = (cf>( M ) + (A25) 

TUB »> ' 

cf (M) - cf (Ai) = - ^ (cf > ( M ) + 2Z#> („)) + cf ( M ) , ( A26) 

where Cq^ (/i) is the Wilson coefficient appearing in the matching of the scalar current in QCD onto HQET 

sb=C^{fi)qh v + --- . (A27) 

Another application of the equations of motion for the vector current id^{q r y^b) — (nib — ni q )(qb) determines this 
Wilson coefficient in terms of those of the vector current as 

C<*V)+^V) = — ^V)- (A28) 

TUB 

At the order we work, the B meson mass can be replaced with the b quark pole mass, and the corresponding mass 
ratios in Eqs. (|A25|) . I|A26|) and (|A27jl are given by 

!^M^ 1 + ^l(_ 61og Ji_4). (A29) 
niB 47r nib 

Combining these relations we find predictions for the Wilson coefficients (/z) , which are confirmed also by explicit 
computation at one-loop order 

DP 0*) = S^L U log JL + 2) , D W (/i) = U log JL + 2 ) . (A30) 

47r \ mb ) 47r \ nib J 

The constraint Eq. i|A25(l can be used to relate the scaling of the 1 + 2D^ (u) / (//) factor to known anomalous 
dimensions. It satisfies the RG equation 

with anomalous dimension 7u(a s ) = — 7t(a s ) — j rn (a s ). We denoted here with jt the anomalous dimension of the 
tensor current defined as 

M^ff+(<? 2 ) = -lt(a s )g+(q 2 ) , lt(a s ) = ^ H (A32) 
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and 7 m is the mass anomalous dimension 

fi—m(fi) = ~/ m (a s )m(fi) , j m (a s ) = -— H . (A33) 

an 7r 

Similar relations among the tensor and axial form factors are obtained starting with the operator identity (valid in 
the NDR anti-commuting 75 scheme) 

id h '(qia ^756) = (m b - m g )<j7 P 7 5 6 - 2qiD ^ 5 b + id^q^b) . (A34) 

Taking the B — > V matrix element gives three relations 

(m| - m 2 v )g+(q 2 ) + q 2 g-(<l 2 ) = (m b - m q )f(q 2 ) - 2d x {q 2 ) (A35) 
~ 9+ (q 2 ) + q 2 h(q 2 ) = (m b - m q )a+(q 2 ) - 2d+(q 2 ) (A36) 
-g-{q 2 ) - {m% - m 2 v )h{q 2 ) = (m b - m q )a_{q 2 ) - 2d_{q 2 ) + s(q 2 ) (A37) 

After using here the l/m b expansions for the c?i ;+ ._ (<? 2 ) form factors, we find the final form of the symmetry relations 
to subleading order in 1 / m b 

(m| - m 2 v )g + (q 2 ) + q 2 g_(q 2 ) = (1 + 2D^^)/C^ ) (ji))(m b - m q )f(q 2 ) - 2df\q 2 ) + ■■■ (A38) 
-g+(q 2 ) + q 2 h(q 2 ) = (1 + 2D^ v \fj,)/C^\fj,))(m b - m q )a + (q 2 ) - 2df(q 2 ) + ■■■ (A39) 
~g„(q 2 ) - {m 2 B - m 2 v )h{q 2 ) = (1 + 2D^ (p)/C^ ( M ))(m 6 - m q )a_(q 2 ) - 2d { °\q 2 ) + ■■■ (A40) 

Together with Eq. I|A23J) . these relations are of phenomenological significance and are used in Sec. V to express the 
contribution of the electromagnetic penguin Q7 to the B — > K*l + i~ amplitude. 

We illustrate in the following the application of Eq. (|A35|I to give an alternative derivation of the power correction 
to a heavy quark symmetry relation presented in |6|]. Consider the combination of form factors 

T(q 2 ) = {m B +m v )g + {q 2 ) + (m B -m v )g-(q 2 ) . (A41) 

The relation Eq. 

gives a prediction for T(q 2 ) at the zero recoil point q 2 nax = my) 2 
,2 ^ ( -\ 1 my - A-m q \ 2 2 (0) 2 



HlL^) = 1 + — 1 HCn) dr (qU) . (A42) 

V m B J m B 

The leading term on the right-hand side was obtained in 0, and the 1 jm b correction was given in <Q] (we correct 
here the sign of the 0(l/m b ) term in the brackets). 
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